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Abstract 

Two-dimensional graphite sheets with a certain type of edges are known to support boundary states localized 
near the edges. Forming a flat band with a sharp peak in the density of states at the Fermi energy, they can 
trigger a magnetic instability or a distortion of the lattice in the presence of electron-electron or electron-phonon 
interactions. We shall discuss a relationship between chiral symmetry, which is the origin of the zero-energy edge 
states, and several types of induced orders such as spin density waves or lattice distortions. We also investigate 
electron correlation effects on the edge states for a wrapped quasi one-dimensional geometry, i.e., carbon nanotube, 
by means of the renormalization group and the open boundary bosonization. 
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1. Introduction 

A fundamental aspect of two-dimensional (2D) 
graphite sheet is its intriguing band structure caused 
by the honeycomb network formed by (T-bonds {sp^) 
on which tt electrons are performing quantum me- 
chanical hoppings. The nontriviality of the band can 
be quantitatively (and mathematically) characterized 
in terms of the Berry (Zak) phase, which is a phase 
degree of freedom in the quantum mechanical systems. 
The importance of the underlying gauge structure 
(the Berry phase) has also been recognized in the con- 
text of the integer quantum Hall effect or a theory of 
macroscopic polarizations[l]. 

A direct consequence of the nontriviality of the band 
is the appearance of edge states which are localized 
near the edges when the system is truncated with a 
certain type of edges. Indeed, Fujita et al. discovered 
peculiar edge states forming a flat band at zero-energy 
(the Fermi energy) in graphite sheets with zigzag or 
bearded (Klein) edges. [2] 
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In the present paper, we consider edge states physics 
caused by electron correlation effects or coupling with 
lattice distortions on the honeycomb lattice with sev- 
eral types of edges. One of our main messages is that 
patters of order parameters [spin density waves (SDW) , 
lattice distortions, etc.] which can be induced in the 
presence of interactions are dictated by a symmetry 
called chiral symmetry that is responsible for the ex- 
istence of the zero-energy edge states (ZES's). Our ar- 
gument is generic and can be applicable to other sit- 
uations such as coexistence of time-reversal symmetry 
breaking order parameters at (110) surfaces of d^2_y2- 
wave superconductors. [3,4] 

As an example, we shall discuss coupling between 
TT electrons and lattice distortions in terms of the Su- 
Schrieffer-Heeger (SSH) model. A possibility of a SDW 
order within the mean fleld theory (MFT) for the Hub- 
bard model defined on the honeycomb lattice with 
edges is also considered. We also give detailed discus- 
sions on electron correlation effects when a graphite 
sheet are wrapped into quasi-lD geometry, i.e., carbon 
nanotube (CNT), where correlation effects are much 
more pronounced due to its low-dimensionality. The 
emergence of magnetic moments at edges is attractive 
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since it opens up a possibility to have a new magnetic 
material made of exclusively light elements, without d- 
or /- electrons. [5] 



2. Origin of zero-energy edge states and a 
Jahn- Teller-like argument 

Typical examples of a truncated graphite sheet 
that supports edge states are those terminated with 
zigzag edges in which edge states appear for +27r/3 < 
Ifcyflol < +1", where ky is the wave number along the 
edges and ao the lattice constant (Fig. 1). On the other 
hand, a graphite sheet with armchair edges does not 
support edge states. The band structure for the zigzag 
case is nontrivial as its Zak phase is equal to tt for ky 
specified above. In the present case, nontriviality is 
also conveniently characterized in terms of a topologi- 
cal object, a set of loops in a parameter space, from 
which we can infer the existence of edges states. [3] 

Although nontriviality of the band implies the ex- 
istence of edge modes, it does not assure that edges 
states are located at zero-energy. Chiral symmetry, 
which means there is no matrix element connecting two 
sublattices in the present case, is another ingredient 
that dictates edge modes to be at zero-energy. Similar 
ZES's arc also known for other systems such as surface 
states of d^2_j,2-wave superconductors, and a states lo- 
calized at a soliton in polyacetylene. All such examples 
can be understood in terms of the topology of bands 
and chiral symmetry in a unified fashion. [3] 

As the edge states form a fiat band (Fig. 1), giving 
rise to a sharp peak in the density of states (DOS) at 
the Fermi energy, a natural question is then what can 
happen in the presence of electron-electron or electron- 
phonon interactions. Once wc understand the origin of 
ZES's, it is possible to have some insights as to what 
kind of orders can occur as follows. Within the mean 
field picture, the band structure near the edges can be 
effectively modified when there arc interactions to lift 
the degeneracy and hence lower the ground state en- 
ergy. However, as ZES's axe "protected" by chiral sym- 
metry, such modifications should be accompanied with 
the breaking of chiral symmetry near the boundaries. 
This is in analogy with the Jahn- Teller theorem in 
which lattice distortions should break crystalline sym- 
metries to lift degeneracy in electronic states. [4] 



3. Coupling with lattice distortions 

As a candidate of interactions that might lift the 

degeneracy of ZES's, let us consider coupling between 
the TT electrons and lattice distortions. We use the fol- 




Fig. 1. (a) A {N, —N) CNT with a zigzag edge. Unit cells are 
shown by broken lines, (b) The corresponding energy spectrum. 



lowing simple model which is an extension of the SSH 
model to 2D graphite [6] : 



HsSH = ( - 1 + Vij) {claCja + h.C.) + Vij, 

(ij)<7 (ij) 

(1) 



where cl^/cia is an electron creation/auiiiliilation op- 
erator at a site i with spin a. Here we treat cr-bonds 
as a classical spring with the tension K and yij rep- 
resents the modulation of the length of a bond with a 
constraint Y2{ij) Vij = 0- We numerically determined 
{Uij } that minimizes the ground state energy, assuming 
the periodicity in ^/-direction with the unit cell shown 
in Fig. (1). Super structures with a larger unit cell can 
be energetically more favorable [7,8], but the present 
choice is enough to demonstrate the relationship be- 
tween chiral symmetry and ZES's. 

When there are no boundaries and hence no edge 
modes, a first order transition with lattice distortions 
occurs at K = Kc ~ 1.4 as we decrease K. (Since the 
DOS is vanishing at the Fermi level, the transition oc- 
curs at a finite K.) When there are zigzag edges, one 
naively expects that the degeneracy between ZES's en- 
hances Peierls instability and a lattice distortion near 
edges occurs at some critical value of is' = K^'^^'^ (< 
Kc). However, this is not possible as chiral symmetry 
is respected for any configuration of {t/ij } in the above 
SSH model. Indeed, ZES's arc stable with respect to 
lattice distortions as shown in Fig. (2). 

For K < Kc, the bulk band structure is completely 
reconstructed by lattice distortions. The corresponding 
loops in a parameter space are trivial and there is no 
ZES even when one truncates the system by zigzag 
edges [Fig. (2)]. 
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4. Correlation efFects and spin polarization 

4.1. Mean field theory for 2D geometry 

Another source of an instability, when there is a flat 
band formed by ZES's, is electron correlation, which 
can give rise to several kinds of orders, a typical ex- 
ample of which is the spin density wave. To model the 
correlation efFects on edge states, let us consider the 
Hubbard model defined on the honeycomb lattice with 
zigzag edges [2] 

-H" = - ^ (c-aCjV + h.c.) rii^rni - /i ^ riia- 



(2) 



We adopt the following mean field ansatz 
HsDw = - ^ (cl^Cja + h.c.) - /i^riio- 

{ij)a icy 

-U ^ {riif} (nil) + f ^ ( (nil) + (nn) riif) (3) 



to investigate possible magnetic structures, where 
(riicr) is dotorminod self-consistently. If there is no 
boundary, the antiferromagnetic order occurs at U = 
Uc (Again, the vanishing DOS at the Fermi energy 
allows the transition to occur at a finite U.). Uc is es- 
timated as 2.23 in MFT whereas Uc ~ 4.5 is obtained 
by the quantum Monte Carlo simulation [9] . 

In contrast to lattice distortions in the previous 
chapter, a finite (ni<r) serves as an on-site potential to 
brealc chiral symmetry and hence lifts the degeneracy 
of ZES's. Indeed, Fujita et al.[2] found that magnetic 
moments are induced near boundaries even when 
U < Uc within MFT. The calculated energy spectrum 
with the self-consistently determined mean field back- 
ground (riic) exhibits the splitting of branches of edge 
modes as shown in Fig. (2). 

4.2. Quasi ID geometry: carbon nanotube 

As we are discussing 2D systems with ID edges, MFT 
is not reliable enough. Also, once we wrap a graphite 
sheet into quasi ID tube geometry, i.e., CNT, the va- 
lidity of MFT is much less obvious due to the quantum 
fluctuation enhanced by the low dimensionality. We 
are now to give detailed analyses for this geometry by 
means of the renormalization group (RG) and the open 
boundary bosonization. [10] Another quasi-lD geom- 
etry, which is a ribbon, was investigated by other au- 
thors. [11,12] Our main focuses are on (i) whether the 
localized charge and spin degree of freedom carried by 
ZES's can escape through the contact between the bulk 
conduction electrons and (ii) whether the total spin at 
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Fig. 2. Top: Band structures of 7r-electrons with lattice dis- 
tortions for (a) if = 1.5 > Kc and (b) AT = 1.2 < Kc- Bot- 
tom: Band structures in the background of SDW mean field 
(a) !7 = 1.5 < Uc and (b) !7 = 2.5 < Uc- 

the boundary is polarized. The flrst question can be 
thought as a variant of the problem of quantum dissi- 
pation, examples of which include the Caldeirar-Legget 
model, the X-ray edge problem, the Kondo effect, etc. 

We consider a thin metallic (N, —N) CNT, where 
there exist one (A'' — 6) or two (A = 9) ZES's when 
truncated with a zigzag edge. Starting from lattice 
models with the Coulomb or the Hubbard interaction, 
we can establish low-energy effective theories which de- 
scribe correlation effects at boundaries. For the case of 
only one edge state, it is given by 

H = Ho + Hi, 

Ho=Hc + €ePe + UeHfUl, 

vpXp 



Hi 



vfX. r r + ^ [s+ J- + j+s-\ , 



(4) 



where He is the Hamiltonian for the conduction elec- 
trons that includes the forward scattering part of the 

Coulomb interaction by the open boundary bosoniza- 
tion (small back scatterings are neglected), J, S repre- 
sent the spin operator of the conduction electrons and 
ZES's, respectively, p,pe the density operator of the 
conduction electrons and ZES's, and n„ denotes the 
number operator for ZES's with spin o. vf represents 
the Fermi velocity and an ultraviolet cutoff Tc is intro- 
duced, which is of order of the inverse bandwidth. Ue 
is the electron correlation between ZES's, ee the chem- 
ical potential for ZES's, and Xp,z,± represents a dimen- 
sionless coupling between the edge states and conduc- 
tion electrons. Bare values of the Kondo-like couplings 
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-^z,± (initial condition for RG analyses) are ferromag- 
netic due to the Fermi statistics; spin polarized states 
do not feel the Coulomb repulsion as the double oc- 
cupancy is automatically prohibited due to the Pauli 
principle, and the kinetic energy is quenched for ZES's 
(Hund's law). 

We perform a perturbative RG analysis up to one- 
loop order by infinitesimally rescaling the ultraviolet 



cutoff, Tc 



From the resultant RG equations. 



it can be seen the bulk conduction electrons and edge 
states are completely decoupled in the infra-red (/ — > 
-|-cxd), Xz^± — > 0. Furthermore, the repulsive bulk inter- 
actions are responsible to suppress the charge fluctua- 
tions at boundaries, and hence doubly occupying a edge 
state becomes unfavorable in the infra-red near half 
filling, e — oo, Ue + 2ee — > -l-co. Then, we conclude 
that edge states do not diffuse into the bulk through 
the coupling with conduction electrons. Our RG anal- 
ysis further indicates that a small deviation of from 
the half-filling (te ~ 0) is irrelevant due to renormal- 
izations by Ap, A^.x- 

In order to see whether edge states tends to form a 
localized moments, we need to consider a thicker metal- 
lic CNT which supports more than two edge states. We 
then consider (9, —9) CNT with zigzag edges for which 
edge states appear for fcyOo = — 87r/9 and +Sn/9. The 
effective Hamiltonian for this case is given by 

Ho = H,+ jpip2 + K,SlS2 + ^ [stS2 + 52+ Sr 
-l-(7e l^nijnix -|-n2tn2ij + €ep&, 



4 -P^f 

+VF\.s'r + ^^^[s+j- + j+s-\, (5) 



where Sa represents the spin operator for ZES's with 
kyao = -|-87r/9 (a = 1) and fc^ao = — 87r/9 (a = 2), 
Pa the density operator, ricra the number operator, and 
S = Si 4- S'2,Pe = pi + P2. Initial conditions are 
given by — K± — —J — —2Ue < 0, te ~ 0, and 
Az = Ax ~ — Ap < 0. RG equations for A^, Ax, 
and hu axe identical to those in the case of one edge 
state. Then, \z and Ax become vanishing in the infra- 
red limit, and charge fluctuations are suppressed e — > 
— CO, Ue + 2ee — > -l-oo. Ou the other hand, RG equa- 
tions for Kz^±^ and 7, which determine the total spin 
carried by the ground state of ZES's, are 



dhi 



2A' 



dl 



h 



A 



dhf 



dl 



A 



(6) 



where Hk^ x '''cKz,± and hi := TcI. We see that the 
Kondo couplings A^.x renormalize the exchange inter- 
actions among edge states A'z,x, making it ferromag- 



netic. That is, the couplings between bulk and edge 
states turn out to assist the spin polarization, leading 
to the ground state with a highest spin S — 1 at the 
boundary. 

The result obtained here is consistent with spin po- 
larization found in the mean-fleld theory in the pre- 
vious section, a density matrix renormalization group 
study for a thin semi-conducting CNT [12] , and an ab 
initio local spin density functional calculation (LSDA) 
[13] for 2D sheet geometry. 



5. Conclusions 

ZES is a hallmark of the nontrivial band of a graphite 
sheet, and is related to rich physics. We have explored 
several types of order that can happen in the pres- 
ence of electron-electron or electron-phonon interac- 
tions, with a main emphasis on chiral symmetry and its 
breaking. The Jahn- Teller-like argument presented in 
the present paper can be potentially applicable to other 
examples. As for spin polarization expected from chiral 
symmetry breaking, we investigated its possibility for 
thin metallic CNT's beyond the mean field treatment. 
Our RG analysis shows that bulk part of interactions 
enhance ferromagnetic interactions among edge states, 
leading to spin polarization at edges. 
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